Space-time Thermodynamics of the Glass Transition 
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We consider the probability distribution for fluctuations in dynamical action and similar quantities 
related to dynamic heterogeneity. We argue that the so-called "glass transition" is a manifestation 
of low action tails in these distributions where the entropy of trajectory space is sub-extensive in 
time. These low action tails are a consequence of dynamic heterogeneity and an indication of phase 
coexistence in trajectory space. The glass transition, where the system falls out of equilibrium, is 
then an order-disorder phenomenon in space-time occurring at a temperature T g which is a weak 
function of measurement time. We illustrate our perspective ideas with facilitated lattice models, 
and note how these ideas apply more generally. 
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A glass transition, where a super-cooled fluid falls out 
of equilibrium, is irreversible and a consequence of ex- 
perimental protocols, such as the time scale over which 
the system is prepared, and the time scale over which 
its properties are observed (for reviews see 0, 0> 13 ) • It 
is thus not a transition in a traditional thermodynamic 
sense. Nevertheless, the phenomenon is relatively pre- 
cipitous, and the thermodynamic conditions at which it 
occurs depend only weakly on preparation and measure- 
ment times. In this paper, we offer an explanation of 
this behavior in terms of a thermodynamics of trajectory 
space. 

Our considerations seem relatively general as they are 
a direct consequence of dynamic heterogeneity [J S 0] 
in glass forming materials. Our primary point is that the 
order-disorder of glassy dynamics is revealed by focusing 
on a probability distribution for an extensive variable 
manifesting dynamic heterogeneity. There are various 
choices for this variable, depending upon the specific sys- 
tem under investigation. We have more to say about this 
later, but to be explicit about what can be learned, most 
of this paper focuses on kinetically constrained models |7j 
and the dynamical actions of these models. 

We show here that due to the emergence of spatial cor- 
relations in the dynamics, i.e. dynamic heterogeneity, ac- 
tion distributions have larger low action tails than what 
would be expected in a homogeneous system. These tails 
indicate a coexistence between space-time regions where 
motion is plentiful and regions where motion is rare. In 
the latter the entropy of trajectories is subextensive in 
time, as previously suggested |8(. The glass transition, 
where the system falls out of equilibrium at long but fi- 
nite observation time, is thus a disorder-order transition 
in space-time. This ordering in trajectory space is not a 
consequence of any underlying static transition (for ther- 
modynamic approaches see 0, E3, El 0)- Further, it 
follows from our results here that the exponential tails 
observed in the aging of soft materials, so-called inter- 
mittency 0,0,03 ( see also ) are a conse- 

quence of dynamic heterogeneity and should be seen in 
mesoscopic measurements in the equilibrium dynamics of 
glass formers. 

Models and distributions of dynamical activity. 



We consider two facilitated models of glass formers to il- 
lustrate our ideas: the lattice model of Fredrickson and 
Andersen (FA) [jjj in spatial dimension d = 1, and its 
dynamically asymmetric variant, the East model [20| . 
The FA and East models serve as a caricatures of strong 
and fragile glass formers, respectively 0,H2 both cases, 
there is an energy function JJ2i n ii where J > sets 
the equilibrium temperature scale, rij is either 1 or 0, in- 
dicating whether lattice site i is excited or not, and the 
sum over i extends over lattice sites. The system moves 
stochastically from one micro-state to another through a 
sequence of single cell moves. In the FA model, the state 
of cell i at time slice t + 1, ra< t+ij can differ from that 
at time slice t, n^t, only if at least one of two nearest 
neighbors, i ± 1, is excited at time t. In the East model 
the condition is that Jij+i,t must be excited. These dy- 
namic constraints affect the metric of motion, confining 
the space-time volume available for trajectories [23|. This 
mimics the effects of complicated intermolecular poten- 
tials in a dense nearly jammed material. Excitations in 
this picture are regions of space-time where molecules 
are unjammed and exhibit mobility. As such, we re- 
fer to rin as the mobility field. For both models, the 
dynamics is time-reversal symmetric and obeys detailed 
balance. The equilibrium concentration of excitations, 
c = (n) = 1/(1 + e J / T ), is the relevant control param- 
eter. The average distance between excitations sets the 
characteristic lengthscale for relaxation, I w c _1 , and 
thus the typical relaxation time, r « c~ a , where A = 3 
for the FA model and A « In c/ In 2 for the East model. 
See Ref. for details. 

The total system under consideration has AT tot sites 
and evolves for t tot time steps, so that the total space- 
time volume is N to t x itot- Within it, we consider a sub- 
system with spatial volume N, and observed time du- 
ration £ Q bs, and use x(t h s ) to denote a trajectory in 
that space-time volume. For the FA and East mod- 
els, this trajectory specifies the mobility field m,t for 
1 < i < N and < t < i bs- The probability density 
functional for x(t h s ), denoted by P[x(t bs)], defines an 
action, £[x(i bs)]; namely, P[x(t ohs )} = exp (-£ [x(t ohs )}). 
In the FA and East models, £{x(t b s )] is a sum over i and 
t of terms coupling n^t to mobility fields at nearby space- 
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FIG. 1: LEFT: Probability distribution of the action P(£) in the FA model at T = 1, subsystem size AT = 60 « 16 t and 
iobs = 320 « 3t. Inset shows probability distribution of number of kinks, PQC) (left and top axes), and of number of excitations, 
P(fC) (bottom and right axes). RIGHT: P(£) but now for t ^ s = 1280 w 12r. The straight line indicates the exponential tail 
exp(£ jii/tobs) (regime a, see text). Inset shows coefficient fi at different T from simulation and theory for both the FA and East 
models. Statistical uncertainties from simulations are smaller than the symbols. 
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in JV x i bs- This extensivity is a general property for 
any system with Markovian dynamics and short-ranged 
interparticle forces. 

The distribution function for the action is 

P{£) ={6(£- S[x(t ohs )})) = n (£) exp (-£) , (1) 

where the pointed brackets indicate average over the 
ensemble of trajectories of length t b s , and fi (£) = 
n (£;N,t b s ) is the number of such trajectories with ac- 
tion £ . When N is much larger than any dynamically 
correlated volume in space, and i bs is much larger than 
any correlated period of time, 



In fi (£, N, Us) = s(£;N, t ohs ) Nt oh 



(2) 



will be extensive in the space-time volume Nt b s . In 
this case, the entropy per space-time point, s(£;N, f bs), 



FIG. 2: Trajectories in the FA model corresponding to the 
regimes (a), (b) and (c) of Fig. 1. Unexcited and excited sites 
are colored white and dark, respectively. Space runs along 
the vertical direction, and time along the horizontal direction. 
Dashed lines delimit the measured space-time volume, iV x 

^obs • 



time points 1 . As such, its average value will be extensive 



1 The microscopic action of either the FA or East models is 



£[x(t)] = £(x t \x t _ St ) + • ■ ■ + £(x St \x a ) - ln[p(xo)]. Here, x t = 
{nu} denotes a configuration at time t, p{xq) is the distribu- 
tion of initial conditions and £{x t j r g t \xt) = — lnp(x t+ g t \xt), 
with p(xt+st\ x t) the probability for an elementary microscopic 
move in time St. In a Monte Carlo simulation, for example, 
St = 1/N and the move is an attempted spin-flip. The ex- 
plicit form of the operator £Jx t+ g t \xt) is given in 0. Here, 
we use continuous time MC l24l where all attempted moves are 
accepted: site i is chosen to flip at time t, jij — ► 1 — m, with 
probability [St) n, + (1 — rii)e~ 1 ^ T ], and time is increased by 
St = (N 1{ + Note- 1 / 7 ")- 1 , where N u and N oi arc the num- 
ber of facilitated up and down sites at time t, respectively. 
The corresponding contribution to the action is £(x t+ g t \xt) = 
— ln<5t + T~ 1 (1 — rii) — In TV (the last term is added to remove 
trivial lniV dependences, also present in standard MC). 
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will be intensive. Similarly, and for the same reasons, 
the mean square fluctuation, %(iV, t a b B ) = (£ 2 ) — (£) 2 , 
is extensive for large enough Nt b s - The onset of this 
extensivity with respect to i bs can be viewed as an order- 
disorder phenomenon, as we will discuss shortly. 

The distribution function for the action can be ob- 
tained from simulation by running trajectories and cre- 
ating a histogram for the logarithm of the probabilities 
for taking steps in the trajectory. Far into the wings of 
the distribution, satisfactory statistics is obtained with 
transition path sampling |2^]. This methodology al- 
lows us to carry out umbrella sampling 26] applied to 
trajectory space 2 . Figure 1 illustrates the P(£)'s we 
have obtained in this way for the FA and East mod- 
els. While the bulk of the distributions is Gaussian, 
for values of the action sufficiently smaller than (£) 
they display exponential tails. Similar distributions are 
found for other measures of dynamic activity, such as 
the number of transitions or kinks in a space-time vol- 
ume, JC = J2i,t \ n i,t (! ~ n Lt+st) + (1 - n itt ) n i<t+st ], or 
the total number of excitations M. = This is 

shown in the inset to Fig. 1 (left) . 

Dynamic heterogeneity and susceptibility. Figure 
2 makes vivid the fact that the exponential tails in these 
distributions are manifestations of dynamic heterogene- 
ity. The extended bubble or stripe in trajectory (a) il- 
lustrates a volume of space-time that is absent of excita- 
tions. Because it extends throughout the pictured time 
frame, the statistical weight for this excitation void is 
dominated by the probability to observe this void at the 
initial time t — 0. As excitations at a given time slice are 
uncorrelated, this probability is Poissonian, exp(— cL), 
where L is the width of the stripe. In the presence of 
a bubble occupying a space-time volume L t bs, the net 
action is £ ss £ b s (N — L) e, where e is the average action 
per unit space-time. In this regime, the probability den- 
sity for the action is then exp [— cL (£)] oc exp (c £ /t Q bs£)- 
This proportionality explains why the slope of the expo- 
nential tail scales linearly with inverse time. The spe- 
cific value of that slope is determined from evaluating 
the action density. Performing the requisite average of 
the action gives e = (£)/Nt b s » —/(c) c lnc, for small 
c, where /(c) rs 2c for the FA model, and /(c) c for 



2 A trajectory of the total system is accepted if the action of the 
subsystem computed from that trajectory lies between £\ and 
£-2- A new trajectory is created from the accepted trajectory by 
implementing shooting and shifting moves l25l . The new trajec- 
tory is accepted if the action of the subsystem again lies between 
£l and £2, and rejected otherwise. The ensemble of trajectories 
created by repeating this process many times provides the distri- 
bution function for the action for the subsystem when that action 
is between £\ and £2 . Yet another part of the distribution is cre- 
ated by performing the sampling again, this time for subsystem 
actions in a different range, say between £3 and £4. Eventually, 
the distribution is mapped out over the entire range of interest 
by splicing together many of these partial distributions. 



the East model. Hence, to the degree that this picture 
is correct, the slope /i/i bs in Fig. 1 should coincide with 
/i = c/e. The inset to Fig. 1 (right) shows that this re- 
lationship holds to a good approximation. Thus, the ex- 
ponential tail in P {£ ) manifests structures in space-time 
that depend upon initial conditions over a time frame 
iobs- It also follows that the entropy for trajectories with 
action that falls in the exponential tails is non-extensive 
in time t b s - 

The tails in P(£) are statistically negligible when i b s 
is large compared to the relaxation time r, but they are 
dominant when t \, s t. The latter regime is where the 
mean square fluctuation or susceptibility, ^obs), is 
non-linear in time. Figure 3(left) shows the growth of 
this quantity with respect to t a bs- The susceptibility per 
unit space-time increases with increasing t b s because in- 
creasing time allows for increased fluctuations. With the 
onset of decorrelation, when i Q b s > t, the susceptibil- 
ity becomes extensive, and x(./V, t ^ s ) / Nt^ becomes a 
constant. This plateau or constant value increases with 
decreasing T because fluctuations become more prevalent 
as c decreases. Indeed, the FA and East models approach 
dynamic criticality as c — > 27]. 

While the plateau value of x (N, t b s ) /Ntobs increases 
with decreasing T or c, the time for the onset to extensiv- 
ity also increases. As a result, if i bs is kept fixed and T is 
varied, \ (N, t Q b s ) will have a maximum. This behavior is 
illustrated in Fig. 3(right). The extremum will approach 
a singularity in the limit of criticality, i b s — > 00, c — > 0. 
For i bs fixed and finite, the extremum is located at a fi- 
nite temperature, T*, namely, r (T*) ~ £ bs, where r (T) 
denotes the equilibrium relaxation time at the tempera- 
ture T. Accordingly, T* is a glass transition temperature. 
In particular, below this temperature, the system can- 
not equilibrate on time scales as short as i bs- In other 
words, with observation times no longer than this, the 
system will have fallen out of equilibrium. For the Ar- 
rhenius FA model, lnr(T) oc 1/T, so that T* in this case 
varies logarithmically with observation time i b s - For the 
super- Arrhenius East model, In r(T) oc 1/T 2 , T* is even 
more weakly dependent upon observation time, going as 
1/vTn £ bs- It is due to this weak dependence on ob- 
servation time that the glass transition is very nearly a 
material property, being confined to a narrow range of 
temperatures. 

Similar behaviors are expected for distri- 
butions of analogous quantities in continu- 
ous force systems. One such quantity is 

2 = Eiev J*o° bB * ex P i ik ■ \ r 3 (t + St ) - T i (*)]}) where 
r*j (t) refers to the position of the jth particle at time 
t, and the sum includes all those particles j that are in 
volume V at time 0. The wave-vector k and time lag St, 
while microscopic, should be such that variation in Q is 
due to diffusion of particles, not simply vibrational mo- 
tion. In the past, studies of dynamic heterogeneity have 
considered the susceptibility for dQ/dt ^ s , regarding the 
time lag St as a variable that can be large [221 l2S] ■ The 
effect of this variability and the differentiation is to focus 
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FIG. 3: LEFT: Action susceptibility per unit length and time, \/Nt obs , as a function of observation time i b s for constant 
temperature T, in the FA model. Inset shows scaled susceptibility, (iJ b 1 s x)/[ um i b 3 -»o°(^ob 1 sX)] as a function of iobs/i~. RIGHT: 
Action susceptibility per unit length and time, x/Nt b s , but now as a function of T at constant t b s - Inset shows free energies of 
trajectories, — In P(£), for observation time t ob s = 3200 at different temperatures. System sizes are N = 16c" 1 . For t ohs = 3200, 
T* ~ 0.44. Statistical uncertainties from simulations are smaller than the symbols. 
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FIG. 4: Free energies of trajectories, Fb(£), for different values 
of the parameter 6. The curve for b = 1 is that of Fig. l(left), 
with T = 1, N = 16c" 1 , t obB = 320. The curves with b = 
1.001, 1.002, 1.003 are obtained by adding (&- 1) to the b = 1 
curve. Triangles show the free energy computed by simulation 
at b — 1.001 using transition path sampling combined with 
free energy perturbation theory. 



a Gaussian random force. In this case, the contribu- 
tion to the action for a transition in infinitesimal time 
dt — > is, up to a constant, — \np[x{t + dt)\x{t)] = 
dt0(t)-f[x(t)}} 2 /AT, which is distributed as -^rf 2 inde- 
pendently of the forces /. The net action for a trajectory 
is simply the accumulation of these infinitesimal steps, 
and will be therefore distributed in a trivial manner. It is 
insensitive to the phenomenon we are after. The action 
for non-linear functions or functionals of the Langevin 
x(t), however, need not be trivial 29]. For variables that 
manifest dynamic heterogeneity, in particular, a mean- 
ingful observable can be constructed out of the probabil- 
ity for finite (not infinitesimal) steps: P[x(t + 5t)\x(t)] = 
J dx t+ dtdx t +2dt ■ ■ ■ p[x(t+dt)\x(t)] ■ ■ ■ p[x(t+St)\x(t+St- 
dt)\, where St is an appropriate coarse-graining time. The 
integrals over intermediate states make the action for fi- 
nite moves, £[x(t + St)\x(t)] = - \nP[x(t + 5t)\x (t)], de- 
pendent on the actual forces /, and the distribution for 
this action will be system dependent. The case of de- 
terministic dynamics is similar in that the microscopic 
action is always vanishing, and a coarse-grained action is 
needed to reveal interesting behaviour. Whether or not 
this coarse-graining is practical, other variables, such as 
Q, are clearly useful and practical. 



on an object that is not extensive in £ ob sj thus obscuring 
order-disorder in space-time. 

Not all variables extensive in space and time are useful. 
This situation is familiar in standard phase transition 
theory. Phase transitions are recognized only through 
physically motivated choices of variables that are then 
shown to reflect the transition. For example, consider 
Langevin dynamics, x(t) — f[x(t)] + ff(t), where rj(t) is 



Order-disorder in space-time. Above T*, the sys- 
tem is equilibrated and disordered with respect to order 
parameters like £ or JC or Q. Below T*, the system is 
ordered with respect to these quantities, and the order 
manifests dependence upon initial conditions. In view of 
the anomalous behavior of \ {N, t hs', T) near T = T* , 
this dependence can be manipulated in a fashion famil- 
iar in the context of equilibrium phase transitions. To 
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this end, it is useful to define b = d In jd£ and a 
corresponding free energy F b (£) = — In 0(f) + b£ = 
— \nP(£) + (b — 1)£. The parameter b is to trajec- 
tory space what reciprocal temperature is to state space. 
In particular, an extremum in \ (T) corresponds with 
a rapid change in {£) with respect to T , which in turn 
implies a bistable free energy, F{,(£), near b = 1 and 

Fig. 4 illustrates this behavior. Coexistence takes place 
at b = b* > 1, the precise value depending upon T and 
i D bs- The free energy for b > 1 can be computed from 
path sampling with a form of thermodynamic perturba- 
tion theory |2g, running trajectories as before (i.e., with 
b = 1) and averaging the functional exp (— A6£[x(i bs)])- 
In the case illustrated, b* = 1.003, very close to the phys- 
ical value 6 = 1. The basin at low action coincides with 
the exponential tail and thus correlation with initial con- 
ditions throughout the observed time frame. 



For an equilibrium system, 6=1. In this case, there 
will be a true dynamical singularity for t D b s — > oo only 
at T = 0. In particular, from the exponential tails of 
P(£), we know that b* = l + /z/t b s , so that b* — ► 1 when 
c — * with t b s /T — > finite. An interesting question 
is whether it is possible to define a dynamical protocol, 
perhaps through some form of external driving, which 
would allow one to tune b from b = 1 to b = 6*, and thus 
observe this dynamical transition at T > 0. 



We thank Hans Andersen and Jorge Kurchan for help- 
ful comments on earlier drafts of this paper. This work 
was supported by the NSF, by DOE grant no. DE-FE- 
FG03-87ER13793, by EPSRC grants no. GR/R83712/01 
and GR/S54074/01, and University of Nottingham grant 
no. FEF 3024. 



[1] Ediger, M.D., Angell C.A., & Nagel, S.R., (1996) J. Phys. 

Chem. 100, 13200-13212. 
[2] Angell, C.A., (1995) Science 267, 1924-1935. 
[3] Debenedetti, P.G., & Stillinger, F.H., (2001) Nature 410, 

259-267 . 

[4] Ediger, M.D., (2000) Annu. Rev. Phys. Chem. 51, 99- 
128. 

[5] Glotzer, S.C., (2000) J. Non-Cryst. Solids, 274, 342-355. 
[6] Richert, R., (2002) J. Phys. Condens. Matter 14, R703- 
R738. 

[7] Ritort, F., & Sollich, P., (2003) Adv. Phys. 52, 219-342. 
[8] Garrahan, J. P., & Chandler, D., (2002) Phys. Rev. Lett. 

89, 035704.1-4 . 
[9] Kivelson, D., Kivelson, S.A., Zhao, X., Nussinov, Z., & 
Tarjus, G, (1995) Physica A 219, 27-38. 
[10] Xia, X., & Wolynes, P.G., (2000) Proc. Natl. Acad. Set. 

USA 97, 2990-2994. 
[11] Franz, S., & Parisi, G, (2000) J. Phys.: Condens. Matter 
12, 6335-6342. 

[12] Bouchaud, J. P., & Biroli, G, (2004) J. Chem. Phys 121, 
7347-7354. 

[13] Cipelletti, L., Bissig, H., Trappe, V., Ballesta, P., & 
Mazoyer, S., (2003) J. Phys. Condens. Matter 15, S257- 
S262. 

[14] Chamon, C, Charbonneau, P., Cugliandolo, L.F., Reich- 
man, DR., & Sellitto, M., (2004) J. Chem. Phys. 121, 
10120-10137. 

[15] Ritort, F., (2004) J. Stat. Mech., art. no. P10016. 

[16] Bramwell, S.T., Christensen, K., Fortin, J.-Y., 



Holdsworth, P.C.W., Jensen, H.J., Lise, S., Lopez, J.M., 

Nicodemi, M., Pinton, J.-F., & Sellitto, M., (2000) Phys. 

Rev. Lett. 84, 3744-3747. 
[17] Buisson, L., Bellon, L., & Ciliberto, S., (2003) J. Phys. 

Condens. Matter., S1163-1180. 
[18] Sibani, P. & Jensen, H.J., (2005), Europhys.Lett., 69, 

563-569. 

[19] Fredrickson, G.H., & Andersen, H.C., (1984) Phys. Rev. 

Lett. 53, 1244-1247. 
[20] Jackie, J., & Eisinger, S., (1991) Z. Phys. B 84, 115- . 
[21] Garrahan, JR., & Chandler, D., (2003) Proc. Natl. Acad. 

Set. USA 100, 9710-9714. 
[22] Jung, Y.J., Garrahan, JR., & Chandler, D., (2004) Phys. 

Rev. E 69, 061205.1-7. 
[23] Whitelam, S., & Garrahan, JR., (2004) J. Phys. Chem. 

B 108, 6611-6615. 
[24] Bortz, A.B., Kalos, M.H., & Lebowitz, J.L., (1975) J. 

Comp. Phys. 17, 10-18. 
[25] Bolhuis, P.G, Chandler, D., Dellago, C, & Geissler, 

PL., (2002) Annu. Rev. Phys. Chem. 53, 291-318. 
[26] Frenkel, D., & Smidt, B., (2002) Understanding Molecu- 
lar Simulation (Academic Press, Boston). 
[27] Whitelam, S., Berthier, L., & Garrahan, JR., (2004) 

Phys. Rev. Lett. 92, 185705.1-4. 
[28] Toninelli, C, Wyart, M., Berthier, L., Biroli, G, & 

Bouchard, J.-P., (2005) Phys. Rev. E 71, 041505.1-20. 
[29] See, for example, Biroli, G, & Kurchan, J., (2001) Phys. 

Rev. E64, 016101.1-14. 



